Nota técnica 


ESTIMACIÓN DE PRECIPITACIÓN DIARIA A 
TRAVÉS DE UN SIG CON IMÁGENES DE RADAR 
METEOROLÓGICO 


e Iván Vilchis-Mata e 
Universidad Autónoma del Estado de México 


e Emmanuelle Quentin + 
Universidad Nacional de Loja, Ecuador 


* Khalidou M. Bá + Carlos Díaz-Delgado + 
Universidad Autónoma del Estado de México 


Resumen 


El presente trabajo está centrado en la integración de información, corrección de 
errores y estimación de datos cuantitativos de precipitación diaria con base en 
imágenes de radar meteorológico banda C. Este desarrollo incluye una propuesta 
metodológica y aplicación geomática en un Sistema de Información Geográfica 
(Idrisi). El caso de estudio se realizó con base en la información generada por el 
radar Cerro Catedral, ubicado en la porción central de México, con una resolución 
espacial y temporal base de 832.78 m y 15 min, respectivamente. En su mayoría, 
los coeficientes de determinación alcanzados (r?) entre la precipitación diaria 
estimada con información de radar y la obtenida en estaciones meteorológicas 
automáticas (EMAs) variaron entre 0.50 y 0.90 al interior de un radio de influencia 
de radar de 150 km. Con base en los resultados obtenidos, la herramienta geomáti- 
ca construida es una alternativa de estimación de precipitación diaria que mejora 
y complementa a la red pluviométrica terrestre. Con ello se facilita la comprensión 
de la variabilidad espacio-temporal de las precipitaciones y justifica su empleo en 
estudios hidrológicos y de gestión integrada de recursos hídricos. 


Palabras clave: precipitación, radar meteorológico, SIG-Idrisi, estimación de 
precipitación. 


Introducción desde los años cuarenta, en el pronóstico me- 


teorológico y en un nivel cualitativo. Hoy en día, 


La evaluación y gestión de recursos hídricos 
en una cuenca o región es una función del 
procesos 
hidrológicos que se llevan a cabo dentro de la 


conocimiento y análisis de los 


misma, y dondelaprecipitaciónesunacomponente 
principal. De manera tradicional, las redes de 
pluviómetros y pluviógrafos han sido utilizadas 
para analizar la variabilidad espacio-temporal de 
la precipitación, y generar un muestreo discreto 
en el espacio y continuo en el tiempo. Por su 
parte, el radar meteorológico ha sido utilizado, 


el radar es una herramienta eficaz en hidrología 
que permite la estimación cuantitativa de 
pluviometría y el conocimiento detallado del 
campo de tormenta, con escala temporal de 
5 o 15 minutos, una resolución espacial de 1 a 
2 km?, en un radio de cobertura hasta de 300 
km. Con ello se obtiene un muestreo continuo 
de la precipitación en el espacio y discreto 
en el tiempo para su análisis sobre amplias 
regiones. Su utilidad en la modelación y gestión 
hidrológica para el pronóstico de inundaciones se 
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adapta perfectamente a la utilización de modelos 
hidrológicos distribuidos (Sánchez-Diezma et al., 
2004). Sin embargo, los datos de radar contienen 
errores asociados con la propagación anómala, 
ecos falsos, bloqueo orográfico y altura del 
haz sobre el terreno (Díaz-Vera et al., 2009). 
Adicionalmente, se puede mencionar la falta de 
una estructura y procesos dentro de un Sistema 
de Información Geográfica (SIG), que faciliten el 
tratamiento de datos espaciales y temporales, a 
fin de disponer de estimaciones pluviométricas 
adecuadas, oportunas y confiables. 

El objetivo principal de este trabajo es 
construir una herramienta geomática en el 
SIG-Idrisi, que permita la integración de datos 
de radar y la generación de información con 
calidad para su uso hidrológico en aquellas 
cuencas insertas en el radio de acción óptimo 
del radar o redes de radares. Los resultados de 
este proceso son evaluados por los coeficientes 
de determinación (1?) entre la estimación de 
precipitación con imágenes corregidas, acu- 
muladas a nivel diario, y los datos obtenidos 
en 36 estaciones meteorológicas automáticas 
(EMAs). Este proceso evalúa la robustez 
e incertidumbre de la presente propuesta, 
que pretende complementar la información 
generada por las redes terrestres de monitoreo 
de precipitación. 


Materiales y método 


El proceso sistemático de corrección de 
errores y análisis de la información generada 
por radar es presentado a través de diez 
etapas y con base en el análisis de imágenes 
PPI (Plan Position Indicator) de resolución 
espacial de 832.78 m, resolución temporal de 
15 min, tomadas a un ángulo de 0.5 sobre 
la horizontal, con un radio de cobertura de 
299.8 km (figura 1) del radar meteorológico 
Cerro Catedral banda C (1= 5.33 cm) ubicado 
en el Estado de México, a los 19* 32* 57” de 
latitud norte, 99* 31 13” de longitud oeste y 
a una altitud de 3 790 msnm (terreno 3 770 
msnm). 


Etapa 1. Obtención de imágenes 


Corresponde a la transferencia de archivos 
crudos PPI del radar Cerro Catedral por parte 
del Servicio Meteorológico Nacional (SMN). 
Su almacenamiento considera por día 96 
imágenes (una imagen cada 15 minutos). Se 
recopiló un total de 17 588 imágenes, para el 
periodo comprendido entre las 00:01 UTC del 
1” de junio de 2008 y las 23:46 UTC del 1? de 
junio de 2009. Es necesario subrayar que no 
todos los días se obtuvieron imágenes cada 15 
minutos; el periodo con más datos ausentes 
se presentó de septiembre a noviembre. 


Etapa 2. Integración de imágenes 


Este proceso se realiza mediante el módulo 
GENERICRASTER (ldrisi). Para ello es nece- 
sario indicar el tipo de archivo (formato de 
1 banda BSQ-band-sequential), número de 
columnas y renglones del formato raster (720 
x 720), coordenadas extremas del radio de 
influencia del radar (299.8 m), así como el 
sistema de referencia geográfico, el cual debe 
ser definido considerando que las líneas de 
azimut (haz del radar) son rectas y los ángulos 
pueden ser usados para representar líneas de 
distancia o rango constante, donde el centro 
de la imagen corresponde con la localización 
del radar. Dado que las imágenes del radar 
están en coordenadas polares, es necesario 
definir y adecuar el sistema de referencia 
USGS Lambert Azimuthal Equal Area (major 
s-ax: 6371000, minor s-ax: 6371000, origin 
long: -99.520280 y origin lat: 19.549300). 


Etapa 3. Giro de imágenes 


La información original presenta un giro 
de 180% debido a la dirección de escaneo 
del radar y es necesario invertir el orden de 
los renglones de la imagen con el módulo 
TRANSPOSE (arisi). 
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Figura 1. a) Ubicación y producto PPI (en formato GIF) generado por el radar de Cerro Catedral y b) EMA's. 


Etapa 4. Conversión de valores de 
reflectividad 


Para la obtención de los valores de dBZ 
(decibelios de Z), se efectúa una conversión 
de los valores dados en el rango de 0-255 (8 
bits) por las imágenes del radar, con el modelo 
propuesto por Environment Canada (Rashedul, 
2005): dBZ = 0.5 * N — 32. Donde N corres- 
ponde al valor dado en el rango de 0-255, y los 
valores resultantes corresponden al rango de 
-32 a 95.5 dBZ. 


Etapa 5. Corrección de errores de estabilidad 


Se refiere a la reducción de efectos de ate- 
nuación del haz emitido por el radar de- 
bido a lluvia sobre el radomo y se efectúa 
con el método desarrollado por Sempere- 
Torres et al. (2002) y Sempere-Torres (2003), 
proporcionando un factor de corrección 
de +3 dBZ de reflectividad cuando existe 
precipitación sobre el radomo y un umbral de 
precipitación que puede afectar la señal igual 
o superior a 7 dBZ (0.1 mm/60 min). 
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Etapa 6. Caracterización de la precipitación 


Corresponde a la identificación de los campos 
de tormenta, grupos de pixeles en valores de 
dBZ de la imagen con base en el algoritmo 
propuesto por Steiner et al. (1995), para zonas 
convectivas (> 40 dBZ) y zonas estratiformes (< 
40 dB2). 


Etapa 7. Conversión en valores de 
reflectividad 


Una vez caracterizada la precipitación, es 
necesario realizar la conversión de dBZ a 


valores de reflectividad (Z) con el modelo Z = 
10 (aBZ/10), 


Etapa 8. Relación reflectividad-intensidad de 
precipitación 


La relación Z-R constituye un proceso 
que facilita la obtención de intensidad de 
precipitación R (en mm/60 min) a partir de la 
reflectividad(Z), paraelloserecurrearelaciones 
empíricas de la forma Z = aR”. En el caso en 
estudio, se utilizó para zonas estratiformes el 
modelo Marshall-Palmer (Sánchez-Diezma, 
2001): Z =200R** (R =(Z /200)%) y para zonas 
convectivas el modelo NEXRAD: Z = 300R'* 
(R = (Z/300)"”*%). Operacionalmente, el radar 
sólo considera el modelo Marshall-Palmer, 
pero como mejora metodológica se utiliza el 
modelo Z-R, en función de la discriminación 
efectuada en la etapa 6. 


Etapa 9. Generación de productos acumulados 


Esta etapa se centra en la estimación de pro- 
ductos cartográficos acumulados de preci- 
pitación, basados en la adición de imágenes 
consecutivas (cada 15 minutos) para un cierto 
periodo de tiempo (horario 4, diario 96 y 
suma de productos diarios para la estimación 
mensual). 


Etapa 10. Corrección por errores orográficos 


Alos productos acumuladosesnecesario corre- 
girles los errores producidos por intercepción 


del haz con el terreno. Para ello, es necesario 
crear una máscara raster, derivada de la 
adición de imágenes consecutivas, a fin de 
identificar los ecos del suelo, zonas situadas 
detrás de un eco de suelo y donde la lluvia es 
minimizada a causa de la energía perdida en la 
intercepción (zonas apantalladas), y los lóbulos 
secundarios causados por la intercepción con 
el terreno de partes de la energía no incluidas 
en el haz principal del radar, que en general 
resultan con un valor “cero”. Adicionalmente 
se elimina el ruido residual por medio de 
un proceso de interpolación y con ello se 
garantiza la estimación de precipitación en 
zonas interferidas por sistemas montañosos 
y edificios volcánicos. En el caso de estudio, 
se identificaron y corrigieron bloqueos oro- 
gráficos del haz, correspondientes a los 
volcanes Jocotitlán, Nevado de Toluca, 
Popocatépetl e Iztaccíhuatl, así como también 
el sistema de laderas NW del volcán Nevado 
de Toluca, y del Ajusco, los cuales generan 
zonas apantalladas en su parte posterior y 
presencia de ecos del suelo debido a lóbulos 
secundarios ocasionando ausencia de datos. 
Finalmente, se realiza un proceso de análisis 
entre los datos pluviométricos estimados con 
el radar (x) y los observados en las EMASs (y) 
con base en el coeficiente de determinación 
(12), es decir, el cuadrado del coeficiente de 
correlación del momento de Pearson (r) de 


2-9 
Dr) "2 (y -y) 


último proceso interpretado como validación 


la forma: = Este 


puntual y no espacio-temporal del método 
aquí propuesto debido a la diferencia de mues- 
treos correspondientes en cada uno de los 
procedimientos de monitoreo pluviométrico. 


Resultados 


El esquema conceptual propuesto, así como 
los procesos descritos, se pueden estructurar 
bajo un esquema geomático, denominado 


Macromodelo en  IÍdrisi, identificando la 
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entrada de datos, el análisis espacial, los 
procesos aplicados a la información y las 
salidas de estimaciones obtenidas. Como 
proceso de validación de la metodología se 
ha considerado como aceptable un valor de 
1? > 0.5. Las estimaciones de precipitación con 
imágenes de radar resultaron adecuadas para 
distancias al radar inferiores a 150 km, con un 


mejor ajuste en aquellos casos donde se puede 


considerar que la precipitación presentó una 
estructura más uniforme (cuadro 1). Fuera de 
este radio, aun cuando se obtuvieron algunas 
estimaciones aceptables, no se observa un 
patrón que garantice las estimaciones. Es 
importante señalar que la diferencia entre los 
valores estimados y el dato puntual observado 
por la red de automática es de + 20-30%, lo cual 
puede ser considerado como una aproximación 


Cuadro 1. Comportamiento del coeficiente 1? a nivel diario para los meses analizados. 


EMA E a a ñ Observaciones 
(julio 08) | (agosto 08) | (abril 09) | (mayo 09) 
1)DF04 0.93 0.69 0.50 0.73 Muy AERIBCRNO, buena explicación del comportamiento de la 
precipitación. 
1)DF05 030 0.21 0.01 0.56 geo o variabilidad en sub y sobreestimación de 
precipitación. 
1)DF06 0.59 0.58 091 0.77 Muy rd buena explicación del comportamiento de la 
precipitación. 
1)DF07 0.53 0.26 0.52 0.86 Significativo, subestimación de precipitación (estratiforme). 
1)GR04 0.28 0.38 — — Poco significativo, EMA inconsistente en suministro de datos. 
1)H101 0.96 0.75 091 0.97 A o buena explicación del comportamiento de la 
precipitación. 
1)HI104 0.76 0.78 ---- 0.23 Significativo, sobreestimación de precipitación (convectiva). 
1)HI06 0.81 0.85 =--- 0.80 Muy significativo, EMA inconsistente en suministro de datos. 
1) 4107 0.14 0.87 0.76 0.85 Son variabilidad en sub y sobreestimación de 
precipitación. 
1)HI108 0.53 0.46 0.96 0.92 Significativo, subestimación de precipitación (estratiforme). 
1)H109 0.70 0.62 0.94 0.25 Sigmtcauyo, variabilidad en sub y sobreestimación de 
precipitación. 
1)HI110 0.75 0.61 =--- 0.55 Significativo, EMA inconsistente en suministro de datos. 
1)Hn11 0.99 0.59 — 0.82 Significativo, EMA inconsistente en suministro de datos. 
1)MO01 0.20 0.83 0.60 0.49 Significativo, EMA afectada por zona de subdetección. 
1)M002 0.53 0.50 0.31 0.46 Significativo, EMA afectada por zona de subdetección. 
1)MX01 0.82 0.75 0.54 0.46 Significativo, EMA afectada por atenuación sobre el radomo. 
1)MX02 0.93 0.21 0.71 0.72 Significativo, subestimación de precipitación (estratiforme). 3 
1)MX03 0.75 0.77 0.95 0.91 Muy significativo, EMA afectada por zona de subdetección. 3 
1)MX04 0.09 0.01 0.14 0.38 Poco SBnlicanyo, EMA afectada por bloqueo orográfico y E 
subdetección. 5 
1)MX09 0.00 0.70 0.61 0.65 Significativo, sobreestimación de precipitación (convectiva). E 
O 
e se E 
1) 0 0.16 o 0.48 0.12 Poco SIgniaaioO, EMA afectada por bloqueo orográfico y 3 
subdetección. E 
1)0001 0.01 0.13 0.97 (0566 Poco significativo, EMA afectada por zona de subdetección. £ 
2)GR02 0.36 0.09 --- ---- Poco significativo, EMA afectada por zona de subdetección. E 
2)GR07 0.21 0.05 0.01 0.02 Poco significativo, EMA afectada por zona de subdetección. E 
2)GTO1 0.36 0.52 0.45 0.70 es e ón variabilidad en sub y sobreestimación de E 
precipitación. = 
ES 
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Cuadro 1 (continuación). Comportamiento del coeficiente 1? a nivel diario para los meses analizados. 


EMA a E a A Observaciones 
(julio 08) | (agosto 08) | (abril 09) | (mayo 09) 

2) HI02 0.93 0.21 0.97 0.61 Significativo, subestimación de precipitación (estratiforme). 

2)MCO01 0.57 0.78 0.72 0.45 Significativo, subestimación de precipitación y zona de 
subdetección. 

2)MC09 0.63 0.04 0.00 0.01 Poco significativo, EMA afectada por zona de subdetección. 

2)PUO1 0.22 0.00 0.58 0.12 [Poco significativo, EMA afectada por zona de subdetección. 

2)PUO3 0.97 0.33 1.00 0.05 Poco significativo, subestimación de precipitación 
(estratiforme). 

2)PUO04 0.16 0.03 0.25 0.15 Poco significativo, subestimación de precipitación 
(estratiforme). 

2)SLO2 0.77 0.80 o 0.92 Muy significativo, no afectada por bloqueo orográfico y 
subdetección. 

2)TLO1 0.04 0.02 0.66 0.33 Poco significativo, subestimación de precipitación 
(estratiforme). 

2)VRO2 0.00 0.01 0.81 0.50 Poco significativo, subestimación de precipitación 
(estratiforme). 

2) VRO7 -—- ---- -—-—- ---- EMA inconsistente en suministro de datos. 

2)VROS 0.01 0.02 0.00 0.00 Poco significativo, variabilidad en sub y sobreestimación de 
precipitación. 


Radio de distancia al radar: 1) inferior a 150 km, 2) superior a 150 km 


de la incertidumbre al momento de utilizar este imágenes del radar meteorológico es aplicable 
tipo de información. a otros radares banda C. Es importante se- 
ñalar que el nivel de corrección dependerá 
Conclusiones y recomendaciones del análisis espacial y temporal, así como de 
las características propias de cada radar. Por 
El método propuesto de estimación de preci- otro lado, la integración de información de 


pitación con datos del radar meteorológico 
de Cerro Catedral con respecto a las técnicas 
tradicionales puntuales (pluviómetros y 


radares aledaños facilitaría un análisis del 
comportamiento de la precipitación, al abarcar 


j e Ñ : zonas de observación mayores al de un solo 
pluviógrafos) de adquisición de la información : es A cid 
A : . radar; esta información podría ser utilizada 
de precipitación proporciona información La dh 
ad a para la vigilancia y el seguimiento permanente 
cuantitativa corregida y detallada, con una , . ay 
, de fenómenos hidrometeorológicos que afec- 
alta resolución espacial (832.78 m) de los NÓ ] 
d tan al territorio nacional. 
campos de tormenta y representa un ejemplo A : as 
día A : : . Entre las aplicaciones hidrológicas de la 
del análisis más complejo de su tipo debido 
e ¿E información obtenida braya la modelación 
a las características orográficas de la zona. A da a 


Se resalta la identificación y clasificación del distribuida de lluvia-escurrimiento, el análisis 


tipo de precipitación, así como la aplicación de riesgo por inundación, así como el balance 
de una doble relación Z-R, para diferenciar hidrológico a nivel de mesoescala, entre 


los rangos de intensidades y, con ello, con- otros. Particularmente, el uso de información 
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siderar las características y evolución de la obtenida del radar Cerro Catedral ofrece una 
precipitación en la zona de influencia del alternativa para estudios hidrológicos y gestión 
radar. Igualmente, el proceso implementado de recursos hídricos en las cuencas Valle de 
en Idrisi de integración, análisis y corrección México y Lerma en la porción correspondiente 
de errores de los datos obtenidos de las al Estado de México. 
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Abstract 


VILCHIS-MATA, 1, QUENTIN, E., BÁ, K.M. €: DÍAZ-DELGADO, C. Daily rainfall 
estimation using a GIS with weather radar imagery. Water Technology and Sciences (in 
Spanish). Vol. IL No. 4, October-December, 2011, pp. 167-174. 


This research focuses on the integration, error correction and quantitative estimation of daily 
rainfall data based on C-Band weather radar imagery. This work includes a methodological 
proposal and a geomatic application in a Geographic Information System (Idrisi). The case 
study was based on the Cerro Catedral weather radar located in central Mexico, with a spatial 
and temporal resolution of 832.78 m and 15 min, respectively. Most determination coefficients 
(1?) for estimated daily rainfall based on radar information and information obtained by 
the Automatic Weather Stations ranged between 0.50 and 0.90, within a 150 km influence 
radius. Based on these results, the geomatic tool built serves as an alternative for estimating 
daily rainfall that improves and complements the rain gauge network. This facilitates the 
understanding of spatial-temporary variability in rainfall and justifies its use in hydrological 
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and studies on integrated water resources management studies. 
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